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Abstract 

In this paper we study some iterative methods for simultaneous approxima- 
tion of polynomial zeros. We give new semilocal convergence theorems with error 
bounds for Ehrlich's and Nourein's iterations. Our theorems generalize and im- 
prove recent results of Zheng and Huang [J. Comput. Math. 18 (2000), 113- 
122], Petko vie and Herceg [J. Comput. Appl. Math. 136 (2001), 283-307] and 
Nedic [Novi Sad J. Math. 31 (2001), 103-111]. We also present a new sufficient 
condition for simple zeros of a polynomial. 
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1 Introduction and notations 

Let / be a monic complex polynomial of degree n > 2. A point £ in C n is said to be 
a root-vector of / if its components are exactly the zeros of / with their multiplicities. 
There are a lot of iterations for simultaneous computation of all zeros of / (see the 
monographs of Sendov, Andreev, Kjurkchiev [1], Petkovic, Herceg, Ilic [2] 
and Kyurkchiev [3]). The famous one is Weierstrass' iteration [4] 

(1) z k+1 = z k -W(z k ), k = 0,1,2, ... , 

where the operator W in C n is denned by W(z) = (Wi(z), . . . , W n (z)) with 

"«*>=nfej < i = 1 < 2 ' -•»>■ 



*This paper is published in: C. R. Acad. Bulg. Sci 59 (2006), No 7, 705-712. 
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Another iteration for simultaneous finding all zeros of / is Ehrlich's iteration [5] 

(2) z k+1 = F(z k ), k = 0,1,2,..., 

where the operator F in C n is defined by F(z) = (F 1 (z), . . . , F n (z)) with 

/(*) 



(3) F(z) = Zl - 



3+1 



Werner [6] has proved that the iteration function F can also be written in the form 

Wi(z) 



(4) F(z) =Zi- 



l + EW 3 (z)/(z t - Zj ) 

3+1 



Ehrlich's method (2) with the iteration function F defined by (4) instead of (3) is known 
as Borsch-Supan's method since in such form it was proposed for the first time by 
Borsch- Supan [7]. The following iteration is due to Nourein [8] 

(5) z k+1 = G(z k ), k = 0,1,2, ... , 

where the operator G in C n is defined by G(z) = (Gi(z), . . . ,G n (z)) with 

WAz) 



(6) Gi(z) =Zi- 



^ + Y.W 3 {z)/{z l -z 3 -W l {z))' 

3+1 



Nourein's method is also known as Borsch-Supan method with Weierstrass' correction. 

Since 1996, some authors [9, 10, 11, 12, 13, 14, 15, 16] have obtained semilocal 
convergence theorems for Ehrlich's and Nourein's methods from data at one point. The 
best results on Ehrlich' method are due to Zheng and Huang [14] and Petkovic and 
Herceg [15]. The best results on Nourein's method are due to Zheng and Huang [14] 
and Nedic [16]. 

In this paper, we present new semilocal convergence theorems for Ehrlich's and 
Nourein's iterations which generalize and improve all previous results in this area. We 
also present a new sufficient condition for simple zeros of a polynomial. The main results 
of the paper (Theorems 2.1, 3.1 and 4.1) will be proved elsewhere. 

Throughout the paper the norm ||.|| p (1 < p < oo) in C n is defined as usual, i.e. 

\\ z \\p = (Sr=i \ z i\ p ) l ^ P ■ For a given point z in C n we define 

(7) d(z) = (di(z),...,d n (z)) and S(z) = min{di(z), . . . , d n (z)}, 

where di(z) = mm{\zi — z 3 \ : j ^ i} for i = 1,2, ... , n. For a given point z in C n with 
distinct components we use the notations 

W(z) (\\\(z) W n (z)\ W(z) 



2 



2 Localization of simple polynomial zeros 

The following theorem is an improvement of a result of Zheng [17]. 

Theorem 2.1. Let f be a monic polynomial of degree n > 2. Suppose there exist 
1 < p < oo and a point z° in C n with distinct components such that 

(9) E(z°) < l/2 l / q and <f>(E(z )) < 1. 

where q is defined byl/p + l/q — 1 and is a real function defined on [0, 1/2 1 / 9 ) by 



(10) <f>(x) 



;i - x){l - 2 1 /ix) V (n - l) 1/p (l - 2 1 lix), 
Then the following statements hold true. 

(i) The polynomial f has only simple zeros. 

(ii) The closed disks 

D l = {zeC: \z-(z° -W t (z°))\ <C\Wi(z°)\}, i = 1,2, . . . ,n, 

where C = 9X/(1 - 0\ 2 ), A = <p(E(z )), 9 = 1- 2 l l q E(z®), are mutually disjoint 
and each of them contains exactly one zero of f. 

Note that under the conditions (9) Weierstrass' method (1) is convergent with the 
second order of convergence (see Proinov [18]). 

3 Semilocal convergence of Ehrlich's method 

Theorem 3.1. Let f be a monic polynomial of degree n > 2, 1 < p < oo andl/p+l/q = 
1 . Define the real function 

1 / 1 \ 71—1 

ax I 1 ax 



1 - (a + l)x)(l - (a + b)x) \ n - 1 1 - (a + b)x 

where a = (n — l) 1 ^ and b = 2 x l q . Suppose that z° is an initial point in C n with distinct 
components satisfying 

E(z°) < l/(a + b) and (f)(E(z )) < 1. 
Then the following statements hold true. 
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(i) Ehrlich's iterative sequence (2) is well-defined and convergent to a root-vector £ of 
f. Moreover, if (f)(E(z )) < 1, then the order of convergence is three. 

(ii) For each k > 1 we have the following a priori error estimate 

ok \(3 fc -l)/2 

( 12 ) V-^\\ P <^-^r\\ W ^)h 

where X = <f>{E{z Q )), 6 = ip(E(z )), A k = /i(£(,2 )A (3fc - 1)/2 ) and the real functions 
tp and /i are defined by 

if)(x) = — — - and ji(x) 



ax i — ax 



(iii) For all k > we have the following a posteriori error estimate 

(13) < ~ n ; 11^(^)11 , 

V V || S|| p _ x _ ^ Afc || v 

where \ k = <p(E(z k )), 9 k = ip(E(z k )) and fi k = fi(E(z k ). 

Setting p = oo in Theorem 3.1 we obtain the following corollary. 

Corollary 3.1. Let f be a monic polynomial of degree n > 2 and let C be a real number 
satisfying 

1 , (n-l)C 2 ( 1-nC \ n 1 

14 < C < and )—- ! . — - -. — < 1. 

v ' ~ n + 1 (1 - nC)(\ - (n + 1)C) \1 - (n + l)C J 

Suppose that z° is an initial point in C n satisfying 

W(z° 



(15) 



d(z c 



< C. 



Then Ehrlich's method (2) is convergent to a root-vector £ of f with the third order of 
convergence. Moreover, the error bounds (12) and (13) hold with p = oo. 

Corollary 3.1 improves Theorem 4.1 of Petkovic and Herceg [15]. Note that 
Petkovic and Herceg have proved only linear convergence of Ehrlich's method under the 
stronger condition ||W / (2: )|| oo < C 5(z°) with C satisfying (14) as well as 

, x „ (n - 1)C 2 (1 + (n - 1)C) ( 1-nC Y 1 ' 1 

(16) f>--= \ (l-Jc)(l-(n-l)C) ( l-(n + l)c ) <l ^ 
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(17) g{0) < 



l-(ra-l)C 
2C 



where g(x) 



l + 2x for0<a;<l/2; 
1/(1 - x) for 1/2 < x < 1. 



Corollary 3.1 shows that Petkovic-Herceg's assumptions (16) and (17) can be omitted. 

The following result improves and generalizes Theorem 1 of Zheng and Huang [14] 
as well as previous results [9, 10, 11, 12]. 

Corollary 3.2. Let f be a monic polynomial of degree n>3,l<p<oo and 1/p+l/q = 
1 . Suppose z° is an initial point in C n satisfying 



(18) 



W(z°) 



d{z°) 



< 



2(n- l)V9 + 2 



Then f has only simple zeros and Ehrlich' method (2) converges to a root-vector £ of f 
with the third order of convergence. Moreover, the error estimates (12) and (13) hold. 

Proof. From Theorem 2.1 and the proof of Corollary 2 of [18] we conclude that / has 
only simple zeros. It is easy to compute that 



<KR) 



{a + l)(a + 2-b) 



1 + 



n-1 a + 2 - b 



n-1 



for R = l/(2a + 2), where a, b and are defined as in Theorem 3.1. By Theorem 3.1 it 
suffices to prove that <f>(R) < 1. If a > e — 1, then 4>(R) < e/(a + 1) < 1. If a < e — 1, 
then 0(i?) < (e — l)/(a — b + 2) < (e — l)/2 < 1 which completes the proof. □ 

Remark. In Corollary 3.2 we consider an initial condition of the type 

E(z°) < 1 



A(n - I) 1 /* + B' 

Note that such initial conditions can be obtained for every A > 1. For example, if A > 1, 
then one can take 

Corollary 3.3. Let f be a monic polynomial of degree n > 2. Suppose that z° is an 
initial point in C n satisfying 



d(z°) 



< 



1 



1.5n + 1.8 



T/ien Ehrlich's method (2) is convergent to a root-vector £ o/ / £/ie £/wrd order o/ 
convergence. Moreover, the error estimates (12) and (13) ZioW with p = oo. 
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Proof. Define <p n = 0(1/(1. 5n + 1.8)), where <fi is defined by (11) with p = oo. The 
sequence (<f) n ) is increasing for 2 < n < 10 and decreasing for n > 10. Hence <p n < 0io < 1 
for n > 2. Now the conclusion follows from Theorem 3.1. □ 

Corollary 3.3 improves Theorem 4.2 of PETKOVIC and HERCEG [15]. Note that the 
authors of this work have proved that Ehrlich's method is convergent under the condition 
|| W(z°) ||oo < C(n)5(z°), where 



C{n) = 



l/(n + 4.5) for n = 3,4; 
1/(1.545 n + 5) for n > 5. 



Corollary 3.4. Let f be a monic polynomial of degree n > 2. Suppose that z° is an 
initial point in C n satisfying 

W(z° 
d(z°) 

where R is the unique positive solution of the equation 



< R = 0.2922.... 

i 



(19) (— ^— J exp— ^— = 1. 

K ' \l-2xj 1 l-2x 

Then f has only simple zeros and Ehrlich' method (2) converges to a root-vector £ of f 
with the third order of convergence and the estimates (12) and (13) hold with p — 1. 

Proof. It follows from Theorem 2.1 and the proof of Corollary 3 of [18] that / has only 
simple zeros. It is easy to show that 4>(x) < g(x) for < x < 1/2, where 4>(x) is defined 
by (11) with p — 1 and g(x) denotes the left-hand side of the equation (19). Therefore, 
4>(R) < 1 which according to Theorem 3.1 completes the proof . □ 

4 Semilocal convergence of Nourein's method 

Theorem 4.1. Let f be a monic polynomial of degree n >2, 1 < p < oo and 1/p+l/q = 
1 . Define the real function 

a 2 x 3 ( 1 a(x — x 2 ) \ n_1 

(1 — mx + bx 2 )(l — (a + 2)x + x 2 ) \ n — 1 l — mx + bx 2 ) 

where a = (n — l) 1 ^, b = 2 1//<? and m — a + b + 1. Suppose that z° is an initial point in 
C n with distinct components satisfying 

E(z°) < 2/(m + Vm 2 -4b) and <j)(E(z )) < 1. 

Then the following statements hold true. 
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(i) Nourein's iterative sequence (5) is well-defined and convergent to a root-vector £ of 
f. Moreover, if (p(E(z )) < 1, then the order of convergence is four. 

(ii) For each k > 1 we have the following a priori error estimate 

ok \ (4 fe -l)/3 

(21) ll^-e|| p <^^A^II^°)l' 

where X = (f)(E(z )), 6 = i/j(E(z )), fi k = /i(£(,2 )A (4fc - 1)/3 ) and the real functions 
tp and /i are defined by 

1 — mx + bx 2 , , . 1 — x 

ip{x) = — -. — — and fi{x) 



1 - (a+ l)x ^ v ' 1 - (a + l)x' 

(iii) For all k > we have the following a posteriori error estimate 

(22) IK -«Mrr£r 11^)11- 

where \ k = <f)(E(z k )), 6 k = ip(E(z k )) and fj, k = fi(E(z k ). 

Setting p = oo in Theorem 4.1 we obtain the following corollary. 

Corollary 4.1. Let f be a monic polynomial of degree n > 2 and let C be a real number 
such that < C < 2/(n + 2 + y/n 2 +4n - 8) and 

r2 - (n - l) 2 C" f l-(n + l)C + C 2 \ n - 1 

1 ' (l-(n + l)C + C 2 )(l-(n + 2)C + 2C 2 ) \l - (n + 2)C + 2C 2 J ' 

Suppose that z° is an initial point in C n satisfying (15). Then Nourein's method (5) is 
convergent to a root-vector £ of f with the fourth order of convergence. Moreover, the 
error bounds (21) and (22) hold with p = oo. 

Corollary 4.1 improves Theorem 2 of NEDIC [16]. Nedic has proved only linear 
convergence of Nourein's method under the condition HW^z^Hoo < C 5(z°) with C 
satisfying 

n n 2 , a 1 - nC 

< C < ; and ft < — — , 

n + 4 + Vn 2 + 8n 1 + (n - 2)C" 

where /3 denotes the left-hand side of (23). 

The following corollary improves and generalizes Theorem 2 of ZHENG and HUANG 
[14] as well as previous results [10, 12, 13]. 
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Corollary 4.2. Under the assumptions of Corollary 3.2 f has only simple zeros and 
Nourein's method (5) converges to a root-vector £ of f with the fourth order of conver- 
gence. Moreover, we have the error estimates (21) and (22). 

Proof. Corollary 3.2 implies that / has only simple zeros. It it is easy to compute that 



cj>(R) 



2a 2 {a + 1) 



(2(a + l) 2 - (2a + l)o)(2a 2 + 2a + 1)) 



1 + 



a(2a+ 1) 



n- l'2(a+ l) 2 - (2a + 1)6 



n-l 



for R = 1/ (2a + 2), where a, b and are defined as in Theorem 4.1. Taking into account 
that b < 2 we get 



(24) 0(i?) < 

and 

(25) 



2a 2 (a + 1) 



(2(a + l) 2 - (2a + l)6)(2a 2 + 2a + 1)) 



1 + 



1 2a + 1 
n — 1 2a 



n-l 



cj>(R)< 



a + 1 



2a 2 + 2a + 1 



1 + 



1 2a + 1 



n — 1 2a 



n-l 



By Theorem 4.1 it suffices to prove that <f>(R) < 1. We shall consider two cases. 
Case 1. Suppose a > 1.8. It follows from (25) that 



(26) 



<I>(R) < g{a) = 



2a 2 + 2a + 1 



exp < i+ £)- 



The function g is decreasing on (0, oo). Therefore, (f>(R) < g(a) < g(l.8) < 1 . 
Case 2. Suppose a < 1.8. From (24) and the obvious inequality b < a we obtain 

2a 2 (a + l) 3 



(27) 



^ < Ha) = (3a + 2)(2a 2 + 2a+l) eXP 2 " 



The function h is increasing on (0, oo). Therefore, 4>(R) < h(a) < h(1.8) < 1 which 
completes the proof. 

□ 

Corollary 4.3. Let f be a monic polynomial of degree n > 2. Suppose that z° is an 
initial point in C n satisfying 



W{z°) 



d(z°) 



< 



1.4n + 2.8 



Then Nourein's method (5) is convergent to a root-vector £ of f with the fourth order of 
convergence. Moreover, the estimates (21) and (22) hold with p = oo. 
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Proof. Define <p n = 0(1/(1. 4n + 2.8)), where is defined by (20) with p = oo. The 
sequence (0 n ) is increasing for 2 < n < 19 and decreasing for n > 19. Hence <p n < 0ig < 1 



for n > 2. Now the conclusion follows from Theorem 4.1. 



□ 



Corollary 4.3 improves Theorem 3 of Nedic [16]. Note that Nedic has proved that 
Nourein's method is convergent under the condition ||rT / (2; )|| oo < 0(11)6(2°) where 



C(n) 



l/(1.64n+ 1.944) for3<n<23; 
1/(1.42 n + 8.7) forn>23. 



Corollary 4.4. Let f be a monic polynomial of degree n > 2. 
initial point in C n satisfying 



Suppose that z° is an 



W(z° 



d(z°) 



<R = 0.2825..., 



where R is the unique solution of the equation 

„3 



(28) 



x 



cxp ■ 



x — X 



[l-3x + x 2 ) 2 ^l-3x + x 



in the interval (0,2/(3 + Vb). Then f has only simple zeros and Nourein's sequence (5) 
converges to a root-vector £ of f with the fourth order of convergence and error bounds 
(21) and (22) with p = 1. 

Proof. By Theorem 3.4 / has only simple zeros. It is easy to show that <fi(x) < g(x) for 
< x < 2/(3 + v^5) , where g(x) denotes the left-hand side of the equation (28) and 
4>(x) is defined by (20) with p — 1. Therefore, 4>(R) < 1 which according to Theorem 4.1 
completes the proof. □ 



References 

[1] Bl. Sendov, A. Andreev, N. Kjurkchiev, Numerical solution of polynomial 
equations, in Handbook of Numerical Analysis, vol. 3 (Eds P. G. Ciarlet, J. L. Lions), 
Amsterdam, Elsevier, 1994, pp. 625-778. 

[2] M. S. PETKOVIC, D. HERCEG, S. Ilic, Point estimation theory and its applications, 
Institute of Mathematics, Novi Sad, 1997. 

[3] N. Kyurkchiev, Initial approximations and root finding methods, Berlin, Wiley, 
1998. 



9 



[4] K. Weierstrass, Neuer Beweis des Satzes, dass jede ganze rationale Funktion einer 
Veranderlichen dargestellt werden kann als ein Produkt aus linnearen Functionen 
derselben Veranderlichen, Ges. Werke 3 (1903), 251-269. 

[5] L. W. Ehrlich, A modified Newton method for polynomials, Comm. ACM 10 
(1967), 107-108. 

[6] W. Werner, On the simultaneous determination of polynomial roots, Lecture Notes 
Math. 953 (1982), 188-202. 

[7] W. Borsch-Supan, Residuenabschatzung fur Polynom-Nullstellen mittels La- 
grange Interpolation, Numer. Math. 14 (1970), 287-296. 

[8] A. W. M. Nourein, An improvement on Noureins method for the simultaneous 
determination of the zeros of a polynomial (an algorithm), J. Comput. Appl. Math. 
3 (1977), 109-110. 

[9] M. S. Petkovic, On initial conditions for the convergence of simultaneous root- 
finding methods, Computing 57 (1996), 163-177. 

[10] M. S. Petkovic, D. Herceg, Point estimation and safe convergence of root- 
finding simultaneous methods, Scientific Review 21-22 (1996), 117-130. 

[11] M. S. Petkovic, S. M. Ilic, Point estimation and the convergence of the Ehrlich- 
Aberth method, Publ. Inst. Math., 62 (1997), 141-149. 

[12] M. S. Petkovic, D. Herceg, Borsch-Supan-like methods: Point estimation and 
parallel implementation, Int. J. Comput. Math. 64 (1997), 327-341. 

[13] M. S. Petkovic, D. Herceg, S. Ilic, Safe convergence of simultaneous methods 
for polynomial zeros, Numer. Algorithms 17 (1998), 313-331. 

[14] S. M. Zheng, Z. D. Huang, On convergence of Nourein iterations for simultane- 
ous finding all zeros of a polynomial. J. Comput. Math. 18 (2000), 113-122. 

[15] M. S. Petkovic, D. Herceg, Point estimation of simultaneous methods for 
solving polynomial equations, J. Comput. Appl. Math. 136 (2001), 283-307. 

[16] J. Nedic, On convergence of Borsch-Supan's method with Weierstrass' corrections, 
Novi Sad J. Math. 31 (2001), No. 1, 103-111. 

[17] S. M. Zheng, On convergence of a parallel algorithm for finding the roots of a 
polynomial, J. Math. Res. Exposition 7 (1987), 657-660 (in Chinese). 



10 



[18] P. D. Proinov, A new semilocal convergence theorem for the Weierstrass method 
from data at one point, C. R. Acad. Bulg. Sci. 59 (2006), No. 2, 131-136. 

Faculty of Mathematics and Informatics 
University of Plovdiv 
Plovdiv 4000, Bulgaria 
E-mail: proinov@pu.acad.bg 



11 



